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^ ■ Abstract 

o ■ 

I Theoretical results for giant resonances in the three doubly magic exotic nuclei ^*Ni, ^°°Sn 

CN . and "2 

Sn are obtained from Hartree-Fock (HF) plus Random Phase Approximation (RPA) 

P5 ' calculations using the DIS parametrization of the Gogny two-body effective interaction. 

^ i Special attention is paid to full consistency between the HF field and the RPA particle-hole 

' ' ' residual interaction. The results for the exotic nuclei, on average, appear similar to those of 

, stable ones, especially for quadrupole and octupole states. More exotic systems have to be 

^ I studied in order to confirm such a trend. The low energy of the monopole resonance in ""^Ni 

^ . suggests that the compression modulus in this neutron rich nucleus is lower than the one of 

' stable ones. 

(N ■ 

O ■ PACS: 21.10.Re, 21.60.Jz, 23.20.Lv 

\o ■ 
o : 

^ ■ 1 Introduction 

O ; 

. Giant multipole resonances (GR) are collective excitations of nuclei that lie at excitation ener- 

gies above the nucleon separation energy (8-10 MeV), have different multipolarities and carry 
different spin-isospin quantum numbers. They have been observed for stable nuclei throughout 
the mass table with large cross sections, close to the maximum allowed by sum rule arguments, 
implying that a large number of nucleons participate in a very collective nuclear motion [1, 2]. 
It is a challenge both to experimentahsts and theorists to study the properties of these states 
for nuclei far from the vahey of stabihty. Not too much has been done from the experimental 
[ side yet: let us just mention the two measurements of the electric dipole GR (GDR) made in 

. neutron-rich oxygen isotopes [3, 4]. Beside GR, there are also low-lying collective excitations, in 

' particular quadrupole and octupole states, which reflect much more than the GR the detail of 

shell structure. More experimental data are available for such states [5] in the case of unstable 
nuclei, giving us information on the modifications of the shell structure far from stability. 

From the theoretical side, more and more calculations of GR and low-lying states are per- 
formed nowadays in the framework of microscopic HF-I-RPA or HFB-I-QRPA approaches. The 
effective nucleon-nucleon interactions used are taken as non-relativistic effective two-body po- 
tentials [6] or relativistic Lagrangians for meson exchange [7]. Such microscopic approaches, 
although less accurate than more phenomenological ones, usually describe reasonably well the 
properties of these states in stable nuclei. 

Among the effective forces used in the non-relativistic approaches, the Gogny force [8, 9] is 
one of those which has been extensively employed for the description of GR and low-lying states 
in doubly closed shell nuclei with the RPA method [10, 11, 12]. Recently, this force has been 
used for the first time in full Quasi-Particle RPA (QRPA) calculations. Chains of isotopes in the 
oxygen, nickel and tin regions have been studied in order to derive the properties of low-lying 
states [13]. 

The purpose of this paper is to present the results of calculations performed in three spherical 
exotic nuclei: ^^Ni, ^'^'^Sn and ^'^■^Sn, and to compare them with those obtained in stable nuclei. 
More precisely, GR and low-lying states in these nuclei will be analyzed and comparisons will 
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be made with systematics and with analogous quantities in the well-known °°Pb. The latter 
nucleus will serve as a reference and, for this reason, results for 208p]-, ^j^j displayed along 
with those of the three exotic nuclei in most Tables and Figures. Let us point out that the results 
presented here for ^°^Pb are new. They have been derived with the DIS parameterization of 
the Gogny force which is the one currently used now. They slightly differ from those of Ref. [10] 
where the older parameterization Dl was employed. 

A point we pay special attention to in the present work is the effect of the full consistency 
of the residual particle-hole (p-h) interaction with the mean field produced by the same force, 
as allowed by the use of consistently combined HF and RPA approaches. In order to analyze 
this effect, we present results where different components of the residual p-h interaction STicli as 
those generated by the spin-orbit or the Coulomb force are switched off. As will be seen, the 
influence of these often omitted components are far from being negligible. 

In the following Section details concerning the parameters of the two-body force, the nu- 
merical methods used for solving the RPA equations are briefly recalled along with a few useful 
formulas. Results are presented and discussed in Section 3. The main conclusions of this work 
are summarized in Section 4. Let us mention that a preliminary account of the present results 
has appeared in the workshop Proceedings of Ref. [14] . 

2 The HF+RPA approach with the Gogny force 

The RPA approach employed here is described in Refs. [11, 10, 12]. The effective force DIS 
proposed by Gogny [8, 9] is used. This finite-range density-dependent interaction describes the 
mean field of the nucleus, and the residual interaction in the RPA calculations is obtained via 
the functional second derivative of the mean field with respect to the one-body density matrix. 
We want to stress that all the terms of the effective force are considered in the HF mean-field and 
in the residual p-h interaction, including the spin-spin component, the Coulomb force and the 
terms produced by the two-body spin-orbit interaction. Only the two-body terms coming from 
the two-body center of mass correction are not included in the RPA matrix elements. Therefore, 
they have been also left out from the mean field calculations. In order to get equivalent binding 
energies and radii, the coefficient of the spin-orbit component of DIS has been reduced from 130 
McV to 115 McV. Such a procedure was previously employed in calculations with the Dl force, 
as explained in Ref. [8]. The Gogny force DIS including this change of the spin-orbit strength 
will be called DIS'. 

In the results presented here, spherical symmetry is imposed. Consequently nuclear states 
can be characterized by their angular momentum J and their parity vr. The individual Hartree- 
Fock wave functions are expanded on finite sets of spherical harmonic oscillator (HO) wave- 
functions containing 15 major shells for all nuclei. For each nucleus, the value of the parameter 
hu of the HO basis is taken as the one minimizing the HF total nuclear energy. 

The RPA equations are solved in matrix form in the p-h representation. RPA energies do 
not appear very sensitive to the value adopted for the HO parameter of the basis. For instance, 
by changing the optimal HF value hoo = 8.7 MeV in ^ospb by 10%, the variation of the ISGMR 
energy (13.46 MeV) is less than .5% and the energy of the first 2+ at 4.609 MeV is changed by 
less than 5 keV. 

Electric transition operators are defined according to: 

A 

QjM = 2X^(1 -Tz{i))jj{qri)YjM{Oi,(t)i), (1) 

i 

where jj is a spherical Bessel function of order J, q a transferred momentum, the third 
component of the nucleon isospin and Yjm the usual spherical harmonics. 

The degree of collectivity of the excited states is measured from their contribution to the 
Energy Weighted Sum Rule (EWSR) 

MiiQjM) = - Eo)\{N\QjM\0)f (2) 

N 

where |0) and \N) are the RPA correlated ground state and excited states, respectively and 
En — Eq their excitation energies. Eq.(2) can also be expressed as the average in the HF ground 
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state \HF) of a double commutator [15]: 



Qjm, H,Qjm \HF). 



(3) 



Therefore, exact values of Mi[Qjm) can be computed from expression (3) whereas smaller 
values will be obtained from (2), reflecting the finiteness of the particle-hole space used in the 
RPA calculations. 

A comparison between the values calculated from (2) and (3) is shown in Figure 1 for '^^Ni 
as an example. As can be seen, with the 15 major shell basis employed, RPA calculations are 
able to describe with a reasonable accuracy the nuclear response for J'^ = 0+, 2"*", 3~, 4+ and 
5~ up to transferred momenta g=1.5 fm~^. 



3 Results 

First, we will discuss the validity of the doubly-magic nature of these exotic nuclei. The single- 
particle neutron spectra obtained in '^^Ni, "'^'^'^Sn and ""^'^^Sn arc shown in Figure 2. The N=50 
gap in '^^Ni and ^'^'^Sn and the N=82 one in ^^^Sn are of the order of 5 McV, which is less than 
20% smaller than the gaps obtained for stable spherical nuclei with same neutron numbers. The 
same is true for the proton gaps at Z=28 in "^^Ni and at Z=50 in tin isotopes. That is, no 
significant reduction of the magic gaps are observed in these nuclei. Therefore, the three exotic 
nuclei are still doubly magic ones and the HF-I-RPA method is applicable to them. 

In what follows, results for states with multipolaritics 0"*", 2"*", 1~ and 3~ are presented for 
four nuclei ^^Ni, and ^'^^Pb, the latter nucleus being included as a reference. 

The strengths shown in the Figures are given in percentage of the EWSR calculated in the 
long wavelength limit q {). The relevant formulas to be used in this limit for the different 
values of J arc given in the appendix of Rcf. [10]. 

In the present calculations the continuum spectrum of the HF Hamiltonian is approximated 
by a discrete one. As a consequence, the RPA strength functions appear in the form of discrete 
peaks. In order to make comparisons with experiments more meaningful, energy centroids will 
be defined in terms of the moments 

Mk [Qjm) = Y.{En - Eo)''\{N\QjMm''. (4) 

N 

of the strength function. Two of these centroids will be used in the following: the mean value of 
the energy Mi/Mq, and the so-called "hydrodynamic" energy -\/Mi/M_i for isoscalar monopole 
resonances. 

As experimental data on GR energies is scarce in exotic nuclei, comparisons will often be 
made with the systematic A~^/^ empirical laws approximately verified in stable nuclei [2]. Values 
from these systematics as well as available experimental data are given in the Tables. 



3.1 Monopole states 

Figure 3 and Table 1 display the results obtained for the Isoscalar Giant Monopole Resonance 
(ISGMR). 

As is well known, the excitation energies of this resonance strongly depends on the compres- 
sion modulus Knm calculated in infinite nuclear matter [16]. One observes in Table 1 that the 

theoretical energies in ^°^Pb, although in good agreement with the empirical 80yl~^/^ law, are 5% 
lower than the experimental value of Ref. [17]. This difference is consistent with the compression 
modulus found in infinite nuclear matter with DIS', i^„m=209 MeV, which is slightly outside 
the interval 220-235 MeV that explains the bulk of experimental data within non-relativistic 
approaches [18]. 

Concerning the three exotic nuclei, we note that resonance energies significantly differ from 
the empirical law only in ^^Ni. It must be noted that, of all three nuclei, ^^Ni is the one where the 

squared neutron-proton asymetry ((A'^ — Z) jAf' most differs from the one of the stable isotope: 
((Af - Z) jAf - {{N - Z) /A)l^^^i^=0.78, 0.36 and -0.23 in '^^Ni, ^^^Sn and lO^Sn, respectively. 
It is therefore tempting to correlate the ~ 1.5 MeV lowering of the ISGMR found in ^^Ni with 
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this large neutron excess, the contribution of the symmetry term K^ym to the finite nucleus 
incomprcssibility Ka being negative [19, 20]. 

The strengths displayed in Figure 3 show that the major part of the EWSR is concentrated 
in a single peak in all four nuclei. This feature explains why the two sets of theoretical energies 
listed in Table 1 are very close to each other. One notes that the fragmentation of the strength 
is almost zero in the N=Z nucleus ^'^''Sn, whereas it is slightly bigger in the other three nuclei 
which have neutron-proton asymmetry {N — Z)/A in the range .21-. 28. 

In Table 2, we show the values of the mean monopole energies Mi/Mq obtained when different 
terms of the residual particle-hole (p-h) interaction are left out of the RPA calculation. Columns 
(1), (2) and (3) refer to the mean energies calculated by leaving out the spin-orbit and the 
Coulomb terms, the Coulomb term and the spin-orbit term, respectively. 

One observes that the spin-orbit part of the residual interaction gives a contribution to 
ISGMR energies ranging from 8% in ''^Ni to 5% in ^"^^Pb. In contrast, the Coulomb contribution 
is larger in Pb (3%) and almost negligible in Ni. These results arc consistent with those discussed 
in Rcf. [18] where '^'^Ca, ^°Zr and ^osp^-, -^^qj-q analyzed with the SLy4 interaction. In the latter 
work, the inclusion in the constrained HF (CHF) of the Coulomb force and of the spin-orbit 
component of the Skyrme interaction was proved to be essential in order to reconcile the value 
of obtained with the Skyrme and Gogny forces. 

3.2 Quadrupole states 

Figure 4 and Tables 3, 4 and 5 display the results obtained for isoscalar quadrupole states. 
Figure 4 shows that in all four nuclei the quadrupole strength is divided essentially between two 
states: the isoscalar Giant Quadrupole Resonance (ISGQR) exhausting ~ 80% of the EWSR 
with an energy in the range 12-16 MeV and a lower-lying state at ~ 3-5 MeV carrying ~ 10%-15% 
of the quadrupole strength. We will label the latter 2'y. 

The theoretical ISGQR energies arc calculated using Mi/Mq excluding the 2f state. The 
results shown in Table 3 are seen to be higher than the A~^/^ systematics by 1.0-1.5 MeV. As the 
latter agrees well with the experimental value in ^°^Pb, it is difficult to draw definite conclusions 
concerning the behaviour of our results in the three exotic nuclei. Let us mention that such large 
ISGQR energies can be understood from a too large spreading of the particle-hole spectrum in 
the 2"*" channel at high energies. Such spreading is a consequence of the value of the effective mass 
of the DIS' interaction (m*/m = 0.7) which is the one giving correct single-particle properties in 
mean-field calculations. As is well known, taking into account the coupling of RPA configurations 
to 2-particle-2-hole (2p-2h) states would reduce this disagreement [21, 22]. Clearly, such a 
coupling should be introduced in the present calculations before reliable predictions for the 
ISGQR in exotic nuclei can be made [23]. Let us mention that the same is true for the other 
giant resonances, with some dependence on the mode quantum numbers [22]. Nevertheless, few 
results have been obtained up to now with such a coupling and it is difficult to foresee the 
magnitude of energy shifts, except for quadrupole and dipole states. 

Our theoretical results for low-lying 2'^ states are presented in Table 4. For these states, 
experimental data exist both for ^"^^Pb [24] and ^'^^Sn [25]. As can be seen, a fair agreement 
between experiment and theory is found in ^'^^Pb and an even better one in -"-^^Sn, with B(E2) 
values being of the same order of magnitude as experimental ones. Let us point out that QRPA 
calculations applied to quadrupole states have been made recently with the DIS interaction for 
a series of tin isotopes including ^^^Sn [13]. In these calculations, the spin-orbit part and the 
coulomb part of the residual interaction were omitted for simplicity reasons. The 2^ energies 
were found larger than the experimental ones by 400 keV in ^''^Sn and 1 Mev in ^^^Sn. The 
corresponding theoretical B(E2) values were lower than experimental ones by at least a factor 
of two. 

These results are consistent with those shown in Table 5 where the same quantities as those of 
Table 4 are displayed. They have been calculated by leaving out from the DIS' p-h interaction 
the spin-orbit and the Coulomb terms, the Coulomb term, the spin-orbit term and no term, 
respectively. One observes that, as previously for monopole vibrations, taking into account the 
spin-orbit part of the residual interaction is essential to get results consistent with experimental 
data. 

Going back to Table 4, 2^ energies are similar in ^^''Sn and ^^^Sn, whereas a comparatively 
low value is predicted in ''^Ni. Let us note that the 2f state in ^^Ni is still higher than the 



4 



one in ^^Ni, the other doubly magic Ni isotope, where the experimental value of the 2^state is 
2.7 MeV and the RPA calculated one is 2.42 McV with DIS'. 

The collectivity of this 2+ state appears larger in ^'^'^Sn than in ^"^^Sn and rather weak in 
^^Ni. Figure 5 displays the transition density ptr of this first 2'^ state in ^^Ni. The definition of 
the transition density is the same as the one given in appendix of Rcf. [10]. One observes that 
the two transition densities arc in phase and that the neutron transition density is higher than 
the proton one and displaced to a larger radius. This mode can therefore be interpreted as an 
isoscalar surface mode dominated by neutron excitation. 



3.3 Dipole states 

Results for the isovector dipole resonance (IVGDR) are presented in Figure 6 and Table 6. 
^'^'^Sn is the nucleus where the giant dipole mode is the least fragmented with 70% of the 
strength concentrated into two peaks. The dipole responses of ^'^'^Pb and ^"^^Sn and to a lesser 
extent of '^^Ni also appear concentrated into two main energy regions. It is expected, that 
the fragmentation is somewhat reduced by the coupling of the RPA modes to 2p-2h states, 
producing smoother strength functions, as in Refs. [26, 27] where Skyrmc forces were used. 

In ^''OSn the mean value Mi/Mq = 19.98 MeV is 3 MeV larger than the systematic 79^"^/^ 
law (17.02 MeV). The EWSR value given in Thomas-Reiche-Kuhn (TRK) unit is 1.59, which 
is large compared to typical experimental values [28]. The IVGDR in ^^^Sn is more fragmented 
than in ^'^''Sn. As in ^'^''Sn the mean energy value, 18.33 MeV, is much larger than systematics 
(79^-1/3 ^ 15 52 MeV) and the EWSR value is 1.58. In the case of ^^Ni, the IVGDR is quite 
fragmented with one major peak and smaller ones at higher energy. The mean energy value, 
20.31 MeV, remains higher than systematics {IQA'^/^ = 18.49 MeV) and the EWSR in TRK 
unit is 1.57. 

It must be said that IVGDR excitation energies calculated with the Gogny force usually 
overestimate experimental data. In the case of ^'^^Pb, the calculated mean value is 16.50 MeV, 
which is quite large compared to experiment (13.43 MeV [28]), but smaller than the result of 
Ref. [10]. Let us note that, ignoring the higher part of the IVGDR response by keeping only 
the strength around the main lower energy peak, considerably improves the agreement with 
systematic estimations : mean energy values become 19.28 MeV, 18.16 MeV, 16.81 MeV and 
14.99 MeV in ^^Ni, ^o^Sn, and ^ospb, respectively. 

In fact, calculated IVGDR energies and EWSR appear quite sensitive to the energy interval 
considered and also to the components of the effective interaction included in the p-h residual 
interaction. This is shown in Table 7 where mean IVGDR energies and EWSR in 208p]-, ^^-^^ listed 
for three energy integration intervals and for RPA calculations where Coulomb and/or spin- 
orbit terms are not included in the RPA matrix elements. One can see that the overestimation 
obtained with the Gogny force decreases by ~ 700 keV when the Coulomb and the spin-orbit 
forces are ignored, which is usually done in RPA calculations employing Skyrme forces, see 
however Ref. 29]. By taking all the terms of the Gogny force and considering the largest 
energy interval, the calculated EWSR given is 1.59 in TRK units. This value is higher than 
the experimental one obtained for a 10-20 MeV energy interval (1-37) [28] but lower than the 
one obtained for a energy interval going up to 140 MeV (1.78) [30]. In this case, however, 
another mechanism, the " quasideuteron effect", is expected to play a major role in the photon 
absorption [30]. 

It is of great interest, beyond nuclear physics itself, to study the amount of excited low-lying 
dipole strength, that is the often called "pygmy" resonances. In terms of EWSR, we obtain much 
less than 1% strength below 10 MeV in Ni and Sn nuclei, and about that amount in ^°*Pb. The 
result for Pb is in agreement with the data of Ref. [31]. The absence of collective states in the 
low-lying region is at variance with the results of relativistic RPA calculations [7] , but agrees with 
the arguing in Ref. [6] . There, it is pointed out that the soft dipole strength should decrease in 
nuclei displaying a neutron skin, compared to that in light halo nuclei because of a more efficient 
coupling to the IVGDR. On the other hand, the coupling to 2p-2h can significantly increase the 
amount of low-lying strength [26, 27]. 

By introducing a very small renormalization factor (1.01-1.03) of the residual interaction the 
isoscalar spurious mode can be made to appear at zero frequency. This factor is introduced 
only in the J'^ = \~ subspace. In Table 8, the values of the energy of this state are shown 
as calculated with or without different parts of the DIS' p-h interaction. For each nucleus the 
same renormalisation factor is used in the four cases. The symbol G 9 means that the RPA 



5 



eigenvalue is imaginary. These results show, as expected, that the consistency between the HF 
field and the residual interaction is important for the treatment of the spurious states. 

3.4 Octupole states 

As shown in Figure 7, the J'" = 3~ states belong to two well-separated energy regions. Only the 
component at energies larger than ~15 MeV can be considered as a genuine giant resonance, the 
High Energy Octupole Resonance (HEOR). Keeping only high energy regions (19-35 MeV for 
^°°Sn, 22-31 MeV for ^32Sn, 22-44 MeV for ^^Ni and 13-28 MeV for ^^^Pb), the mean calculated 
HEOR energies are 28.16 MeV, 26.06 MeV, 29.51 MeV and 23.20 MeV, respectively. These 
values give systematica EqA^'^/^ , with Eq = 130, 132, 126, and 137 in the four nuclei, to be 
compared with the usual estimate 1\QA~^/^ [32]. Previous studies in stable nuclei [10] gave 
values between 130A~^/^ and 140^4"^/^ for heavy nuclei and around 120^~^/^ in lighter ones. 
Wc therefore do not observe a strongly different behaviour of HEOR energies in exotic nuclei 
compared to the one previously obtained along the valley of stability. 

The characteristics of the low energy 3~ states are reported in Table 9. The influence of the 
different components of the DIS' force included in the p-h interaction is also shown. The effect 
of the spin-orbit term appears to be smaller than for the quadrupole states in Table 5, especially 
for ^^Ni. 

3.5 Isovector strength 

In Figures 8-10, the fractions of the isovector EWSR carried by the J'" = 0"*", 2+, 3~ states is 
drawn. In this case, systematics for stable nuclei are not yet well known [2] and is not reported. 
Note that only the transition operator is changed compared to the isoscalar case in Figures 3, 
4 and 7. From the comparison between the two sets of figures, a much larger fragmentation of 
the strength is found in the isovector case, and a mixed (isoscalar-isovector) character of several 
states appears, as expected, in particular in ^^Ni. 

4 Conclusion 

To summarize, we have presented the results obtained for different giant resonances in three 
doubly magic exotic nuclei, using the HF-I-RPA approach and the Gogny force. The largest 
difference with usual doubly magic nuclei inside the valley of stability occurs in ^^Ni where 
the ISGMR appears significantly lower than systematics. This seems to be due to the large 
proton- neutron asymmetry of this nucleus. 

The fragmentation of the isovector dipole strength has to be explored further in order to see 
the correlation or the no-correlation with proton-neutron radius differences. In particular, the 
nature of the double-peaks obtained in tin isotopes remains to be determined. 

Results obtained in the three exotic nuclei for the ISGQR and HEOR resonances are similar 
to those of ^°^Pb, but more exotic systems have to be studied to confirm such a trend. 

Low energy states and B(E2) values appear to be well reproduced within the present ap- 
proach, in particular the first 2^ in ^^^Sn. 

From a more general point of view, we have found that the spin-orbit component of the p-h 
residual interaction plays a very important role in the structure of the low-lying quadrupole and 
octupole states, as it strongly influences both excitation energies and transition probabilities. 
Similarly, our results show that including the Coulomb force in the RPA p-h matrix elements 
significantly affects IVGDR energies and EWSR. 
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Figure 1: Comparisons between the exact EWSR of eq.(3) (solid line) and those deduced from 
eq.(2) (dotted line) in ^^Ni for the RPA states with J'" = 0+, 2+, 3", 4+ and 5". The unit of 
the EWSR scale is e^ MeV. The abscissa q represents the transferred momentum 
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Figure 2: Single particle levels in the vicinity of the Fermi surface for neutrons in the three stud- 
ied exotic nuclei. Filled and empty levels are represented by full and dashed lines, respectively. 
The labels indicate the quantum numbers {nlj) of the levels. 
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Figure 3: Fraction of the EWSR carried by isoscalar J'^ = 0"^ states in the four studied nuclei. 
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Figure 4: Fraction of the EWSR carried by isoscalar J'^ = 2+ states in the four studied nuclei. 
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Figure 5: Neutron (full line) and proton (dashed line) transition densities for the first 2+ state 
in ^^Ni. 
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Figure 6: Fraction of the EWSR carried by isovector J'^ = 1~ states in the four studied nuclei. 
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Figure 7: Fraction of the EWSR carried by isoscalar = 3 states in the four studied nuclei. 
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Figure 8: Fraction of the isovector EWSR carried by J'^ = 0"*" states in the four nuclei. 
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Figure 9: Fraction of the isovector EWSR carried by J"^ = 2+ states in the four nuclei. 
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Figure 10: Fraction of the isovector EWSR carried by J'^ = 3 states in the four nuclei. 



0+ r=o 


Ml 
Mo 


/ Ml 
Vm_i 


80 


Exp 


7«Ni 


17.17 


17.07 


18.72 






17.22 


17.18 


17.23 




i32Sn 


15.29 


15.22 


15.72 




208 p]-, 


13.46 


13.42 


13.50 


14.17± 0.28 



Table 1: Mean values and "hydrodynamic" centroids of ISGMR energies in MeV obtained with 
the DIS' force in the four studied nuclei compared with the empirical SOA"^/^ law and the ^^^Ph 
experimental value from Ref.[17]. 



Ml/Mo 


(1) 


(2) 


(3) 


(tot) 


78Ni 


18.55 


17.10 


18.59 


17.17 




18.19 


16.81 


18.54 


17.22 


132 Sn 


16.07 


15.06 


16.26 


15.29 


208p|g 


13.73 


13.05 


14.10 


13.46 



Table 2: Mean ISGMR energies in MeV obtained by leaving out from the DIS' p-h interaction: 
(1) the spin-orbit and the Coulomb terms, (2) the Coulomb term, (3) the spin-orbit term, (tot) 
no term. 
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ISGQR 


DIS' 


64A-V3 


Exp. 




15.94 


14.98 




lOOSn 


15.13 


13.79 






13.79 


12.57 




208 p]^ 


11.98 


10.80 


10.60 



Table 3: Mean values of ISGQR energies in MeV obtained with DIS' in the four studied nuclei 
compared with the empirical 64^1"^/^ law and the 208p|-, experimental value from Ref. [2]. 



2+ 


E 


B(E2) 


Experiment 
E(MeV) 


B(E2) 


V8Ni 


2.73 


466 






lOOSn 


3.84 


1431 






i32Sn 


3.97 


1134 


4.041 


1400 (600) 


208pb 


4.609 


2781 


4.08 


3180 (160) 



Table 4: Energies in MeV and corresponding B(E2) in e^fm^ of 2^ states calculated with the 
DIS' interaction. Existing experimental data from Refs. [24] and [25] are also listed. 



2t 


(1) 




(2) 




(3) 




(tot) 






E 


B(E2) 


E 


B(E2) 


E 


B(E2) 


E 


B(E2) 


V8Ni 


3.53 


257 


2.84 


456 


3.43 


271 


2.73 


466 


lOOSn 


4.64 


1103 


3.95 


1552 


4.48 


1041 


3.84 


1431 


Sn 


4.61 


775 


4.04 


1182 


4.53 


770 


3.97 


1134 


208pb 


5.15 


2305 


4.65 


3145 


5.09 


2123 


4.61 


2781 



Table 5: Energies in MeV and B(E2) of 2f states obtained by leaving out from the DIS' p-h 
interaction: (1) the spin-orbit and the Coulomb terms, (2) the Coulomb term, (3) the spin-orbit 
term, (tot) no term. 
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IVGDR 


DIS' 


79^-1/3 


Exp. 




20.31 


18.49 




lOOSn 


19.98 


17.02 






18.33 


15.52 




208 p]^ 


16.50 


13.33 


13.43 



Table 6: Mean values of IVGDR energies in MeV obtained with DIS' in the four studied nuclei 
compared with the empirical 79A~^/^ law and the 208p|j experimental value from Ref. [28]. 



2U8p^ 


(1) 

<E> EWSR 


(2) 

<E> EWSR 


(3) 

<E> EWSR 


(tot) 

<E> EWSR 


Exp. 
EWSR 


[0 - 140] 
[0 - 20] 
[10 - 20] 


15.88 1.63 
15.10 1.41 
15.20 1.39 


15.70 1.62 
15.31 1.47 
15.17 1.49 


16.71 1.59 
15.83 1.33 
15.90 1.32 


16.50 1.59 
15.86 1.42 
15.95 1.41 


1.78 
1.37 



Table 7: Mean IVGDR energies in MeV and EWSR in TRK units for ^"^Pb calculated by 
leaving out from the DIS' p-h interaction: (1) the spin-orbit and the Coulomb terms, (2) the 
Coulomb term, (3) the spin-orbit term, (tot) no term. The three lines show the results obtained 
for the three energy intervals given in MeV in the leftmost columm. The rightmost column gives 
experimental EWSR in TRK units. 



I7p T=0 


(1) 


(2) 


(3) 


(tot) 


l^^Sn 






2205.78 


4.26 


2U8pb 


G 9 


G 9 


1605.19 


2.29 



Table 8: Energy in keV of the isoscalar 1^ spurious state calculated by leaving out from the 
DIS' p-h interaction: (1) the spin-orbit and the Coulomb terms, (2) the Coulomb term, (3) the 
spin-orbit term, (tot) no term. The symbol G 9 means that the RPA eigenvalue is imaginary. 



3r 


(1) 

E B(E3) 


(2) 

E B(E3) 


(3) 

E B(E3) 


(tot) 
E B(E3) 


Exp 
E B(E3) 


v«Ni 
loogn 

i32Sn 

208pb 


7.95 0.170 
7.26 0.130 
5.78 0.123 
3.55 0.725 


7.80 0.221 
6.95 0.149 
5.60 0.139 
3.38 0.782 


7.87 0.181 
7.13 0.128 
5.72 0.124 
3.57 0.677 


7.70 0.231 
6.82 0.147 
5.53 0.140 
3.39 0.727 


2.6 0.611 (120) 



Table 9: Energies in MeV of the first 3~ state and corresponding B(E3) in lO^e^/m^ calculated 
by leaving out from the DIS' p-h interaction: (1) the spin-orbit and the Coulomb terms, (2) the 
Coulomb term, (3) the spin-orbit term, (tot) no term. Experimental data from Ref. [33] is also 
listed. 
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